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A set of exact quasi-local conservation equations is derived from the Einstein's equations using 
the first-order Kaluza-Klein formalism of general relativity in the (2,2)-splitting of 4-dimensional 
spacetime. These equations are interpreted as quasi-local energy, momentum, and angular momen- 
tum conservation equations. In the asymptotic region of asymptotically flat spacetimes, it is shown 
that the quasi-local energy and energy-flux integral reduce to the Bondi energy and energy-flux, 
respectively. In spherically symmetric spacetimes, the quasi-local energy becomes the Misner-Sharp 
energy. Moreover, on the event horizon of a general dynamical black hole, the quasi-local energy 
conservation equation coincides with the conservation equation studied by Thorne et al. We discuss 
the remaining quasi-local conservation equations briefly. 

PACS numbers: 04.20.Cv, 04.20.Fy, 04.20.-q, 04.30.-w; quasi-local conservation law, quasi-local 
observable, black hole, Kaluza-Klein theory 



I. INTRODUCTION AND KINEMATICS 

In general relativity there have been many attempts to obtain quasi- local conservation equations ■ One of the 
motivations of these efforts is that the quasi-local conservation equations allow us to predict certain aspects of the 
future of a quasi-local region of a given spacetime without actually solving the Einstein's equations for that region. 
Recall that in the Newtonian theory, the conservation of momentum XI P = constant immediately follows from 
Newton's third law, which is no more than the consistency condition implementing Newton's second law. In general 
relativity, the consistency conditions for evolution are already incorporated into the Einstein's equations through the 
constraint equations, from which global conservation equations were found. The purpose of this letter is to show that, 
from the Einstein's equations, one can find conservation equations of a stronger form, namely, quasi-local conservation 
equations. These equations, interpreted as quasi-local energy, momentum, and angular momentum conservation 
equations, are exact and unique in the sense that they are obtained by integrating the Einstein's "constraint" equations 
over a compact two-surface [Q. 

Let us start from the following line element 

ds^ = -2dudv - 2hdu'^ e^pab {dy" + A^du + AJ'dv) {dy^ + A^du + A^dv) , (1) 

where +, — stands for u^v, respectively. The geometry can be understood as follows. The hypersurface u = 
constant is an out-going null hypersurface, and the hypersurface v = constant is either timelike, null, or spacelike, 
depending on the sign of 2h. The intersection of two hypersurfaces u,v = constant defines a spacelike compact 
two-surface N2, on which we introduce the coordinates y°'{a — 2, 3). The metric on N2 is written as 

<^ab = e'^pah, (2) 
where pab is the conformal two-metric satisfying the condition that 

det Pab = 1, (3) 
and e'^ is the area element. Notice that v is the affine parameter of the null vector field 

which rules the out-going null hypersurface u — constant. If we further assume that j4_° — 0, then the metric 
(|l|) becomes identical to the Newman-Unti metric [Q. In this letter, however, we shall retain the A_^ field, since 
its presence will make the Af2-diffeoniorphism invariant Yang-Mills type gauge theory aspect of this Kaluza-Klein 
formalism transparent. Let us mention that the coordinates used in the above construction arc not unique, which 
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means that there are residual symmetries that preserve the metric (|l|). These residual symmetries consist of the 
diffeomorphisms of N2, the reparametrization of u, and the shift of the origin of the affine parameter v at each point 
of -/V2, and were studied in detail in Jist . 

The spacetime integral /o of the scalar curvature of the metric (0) can be written as 

/q = jdu dv d^y Lq + surface integral, (5) 



where Lq is given by |16|,|11|-[13[ 

Lo - -h^''pabF+^F^^+e''{D+a){D^a) - ^e W'(i5+Pac)(i?-P6rf) + e'^i?2 
-2e''{D^h){D^<j) - hc-{D^af + i/ie W'(^-Pac)(i?-PM). (6) 

Here R2 is the scalar curvature of iV2, and we defined the diffA''2-covariant derivatives as follows, 

= d+A^ - d-A_^ - [A+,A^]l, (7) 
D±a^d±G-[A±,a]i^, (8) 
D±h^d±h-[A±,h]i^, (9) 

D±Pab = d±pab ~ [A±,p]hab, (10) 

where d+ = d/du, 9_ = d/dv, da = d/dy"^, and [A±, *]l is the Lie derivative of * along the vector field A± := A^da- 
In addition to the eight equations of motion that follow directly from the integral Iq by the variational principle, 
there are two supplementary equations associated with the partial gauge-fixing of the general metric to the metric 
(|l|). These equations, from which two quasi- local conservation equations are shown to follow, are obtained by varying 
the Einstein-Hilbert action before the partial gauge-fixing condition is introduced. Here we present them without 
derivation, and they are given by [ p^ 

(i) - e^Dla - ^e^iD+af - {D^h){D+a) + {D+h){D^a) 



+2he^{D^h)iD^a) + c'^F^^Oah - ^e'' p'^" p'=''{D+pac)iD+pbd) + daip"'dbh 



+te^{i?2 - {D+a)iD-a) + ^p'^" p'=''iD+pac){D^pbd) + \e'' PabF+^F^^} 

We'^{{D^af ~ ipV'P-Pac)(i?-P6d)} = 0, (11) 
(ii) e^D+D^a + c^D^D+a + 2g'' {D+(j){D^(j) - 26" {D^h){D^a) 

-^e^-p^,F+-F+b _ g.^^ _ /je-{(i^_a)2 - ipV'(^-Pac)(i?-PM)} = 0. (12) 



II. A SET OF QUASI-LOCAL CONSERVATION EQUATIONS 

Notice that the equations (|ll]) and (|l2|) are first-order in derivatives, so that they may be regarded as two 
"constraint" equations. Thus, in this formalism, the natural vector field that defines the evolution is 9_. Then the 
momenta tt/ = {tt^i, tTo-, tTq, tt"''} conjugate to the configuration variables = {h, a, pab} are defined as 

and are given by 

TTh = -2e"(i^-a), (14) 
-K^ = -2e'"(L>_/i) - 2hc''{D^(j) + e'"(L>+cr), (15) 

TTa = e2"pafcF+^, (16) 

= he^p'''p'"'{D^p,a) - ie>'^V'"'(^+Pcd). (17) 
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Notice that traceless 



due to the identities 



The "Hamiltonian" function Hq defined as 



PabTT"''' - (18) 
p''''D±Pab = 0. (19) 

Ho := iTidW - Lo (20) 



is found to be 

Ho = H - Aj'Ca + surface terms, (21) 

where H and Ca are given by 



+ l7r„(i?+a) + l-rr-^D+pab) + ^e'^ p'^'p^'iD+PacKD+pM) - e'^i?^, (22) 

Cq = d+TTa - dbiA^TTa) - TTbdaA^ - TT^dad + OaTTcr " T^hdah - TT^^daPbc 

+96(7r''Vac) + d,{T:^-pab) - ajTr^Vfec). (23) 
In terms of these canonical variables {tt/, (7^}, the supplementary equations (|ll|) and (|l^) can be written as 

(i') TT^^D+Pab + TTaD+a - hD+TTh - 8+ [h TTh + 2e''D+cj) 

+da {h TThA+ + 2A^c''D+(j + 2he-'' p^^iTb + 2p''''dbh^ = 0, (24) 
(ii') H - d+TTh + da (a/tt,, + e-"p"^b) - 0. (25) 
Moreover, we have two more first-order equations 

Ca = 0, (26) 

which follows trivially by varying Hq with respect to A_^ (or by varying the Einstein-Hilbert action [|l^,^ , which re- 
quires rather lengthy computations). The equations ( |2^), (p5|), and ( ^) are the four Einstein's "constraint" equations 
in the gauge (|l|). Notice that the equations (^) and (p5[) are divergence-type equations. If we contract the equation 
( p6| ) by an arbitrary function i^" such that 

d±C = 0, (27) 
then the resulting equation is also a divergence-type equation, 

n''''£^Pab + Tra£i<7 + T^h£^h + ^a£iA^ - O+i^T^a) 

+da ( - + 27r'^''f + A^S,\b) = 0, (28) 



where £^ is the Lie derivative along ^ := ^°9a. 

The integrals of these equations over a compact two-surface N2 become, after a suitable normalization, 

^U{u, v)^-^ (fd^y U'^^'D+pab + T^aD+a - hD+iTh) , (29) 

^PM^^fd^H, (30) 

^i(u, jd^y (^°'i^«Pab + i^o£^<J - h£^^h - A^£^T,a) , (31) 
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where U{u,v), P{u,v), and L{u,v]^) are defined as 



U{u, y)-^^^ {h T^h + 2e"Z?+a) + U\ (32) 

P{u,v):^^<l^d'y{n^) + P'>, (33) 

L{u,v;0-^^<fd'y{C7Ta) + L''. (34) 

Here H is the Hamiltonian function given by (p^), and C/°, P°, and are undetermined subtraction terms, which 
however must be u-independent, 

du du du 

in order to satisfy the equations (^, (|3^), and (^, respectively. 

One can write the r.h.s. of the equation (^9|) in a more symmetric and suggestive form as follows. Let us contract 
the equation ( p6| ) with and integrate the resulting equation over N2 to obtain the following equation, 

jd^y (A^d+^a) = jd^y (tt'''' £ A^Pab + tt„£a+^ ~ h£A^^h) ■ (36) 

If we use the definitions of diffA^2-covariant derivatives D± and the equation (^), then the equation (29) can be 
written as 

^f/(w, = jd^y {^""d+Pab + ^ad+n - hd+TTh - A^d+TTa) . (37) 

Notice that the r.h.s. of the conservation equations (|3|) and ( ^7|) match exactly, if we interchange the derivatives 

£^ < — > d+ (38) 

in the two-surface integrals. 

In a region of a spacetime where d/du is timelike, these quasi-local conservation equations relate the instantaneous 
rates of changes of two-surface integrals at a given u-time to the associated net flux integrals, and they form a 
"complete" set of quasi-local conservation equations since they follow directly from the four Einstein's "constraint" 
equations. Let us remark that, unlike the Tamburino-Winicour's quasi-local conservation equations |^ which are 
"weak" conservation equations since they assumed in the derivation the Ricci flat conditions (i.e. the full Einstein's 
equations), our quasi- local conservation equations are "strong" conservation equations since we used the four Einstein's 
"constraint" equations only. 



III. QUASI-LOCAL ENERGY AND ENERGY-FLUX INTEGRAL 



The equations (29), (|3C), and ( |3l| ) are in fact quasi- local energy, momentum, and angular momentum conservation 
equations, respectively In this section we shall focus on the quasi- local energy conservation equation (p9|), and 

defer discussions of other conservation equations to the section IV. In order to define a quasi- local energy associated 
with a given two-surface N2, we have to introduce a subtraction term referring to that region only. In general the 
subtraction terms for a quasi-local energy are not unique, and the "right" subtraction term may not even exist at all 
in a generic situation. One natural criterion is that the subtraction term must be chosen such that the quasi-local 
energy reproduces "standard" values in limiting cases. One possible candidate is 



A 

167r' 



(39) 



where A is the area of N2- However, this subtraction term introduces a restriction on the admissible two-surfaces N2 
or, equivalently, on the vector field d/du, due to the condition 



OA 
du 



= 



(40) 
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that follows from (p5[). But this is just the condition that the Bondi time function satisfies at the null infinity. 



Therefore the time function u satisfying the condition (40) may be regarded as a quasi-local, finite distance analog of 



the standard Bondi time function. The condition (^ can be interpreted as follows: In order to evaluate and compare 
quasi-local energies of a given two-surface 7V2 at two different times in a physically meaningful way, there must exist 
certain requirements on N2, and the minimum requirement is that the area of the two-surface N2 is independent of u. 
The quasi-local energy integral for a two surface with the area A is then given by 



v) = ^f d'y {h + 2e^D+a) + ^ A, (41) 



which in general is not positive-definite, and does not possess monotonicity either. The lack of these properties seems 
to be related to the fact that gravitational binding energy is always negative. However, as we shall see shortly, the 
quasi-local energy function at the null infinity reduces to the Bondi mass, which is positive-definite and decreases 
monotonically as the time u increases. 



III-l. The Bondi mass-loss formula 

In asymptotically flat spacetimes, the metric becomes 



ds^ — > -2dudv - ( 1 )du^ + v^dn^ (42) 



V 



as w — > 00. Thus d/du is asymptotic to the timelike Killing vector field. In general, the energy-flux across a 
two-surface is given by the energy- momentum tensor Tq^, which is of the form 

Tor ~ n^d,(j). (43) 



The integrand of the r.h.s. of (37) is of this form, and we therefore expect that it represents the energy-flux carried 
by gravitational radiation crossing N2. Then the l.h.s. of (37) should be the instantaneous rate of change in the 
gravitational energy of the region enclosed by A^2- The energy-flux integral in general does not have a definite sign, 
since it includes the energy-flux carried by the in-coming as well as the out-going gravitational radiation. But in 
the asymptotically flat region, the energy-flux integral turns out to be negative-definite, representing the physical 
situation that there is no in-coming flux coming from the infinity. 

Let us now show that the equation ( p9| ) reduces to the Bondi mass-loss formula [ p^ in the asymptotic region 
of asymptotically flat spacetimes. The asymptotic fall-off rates of the metric coefficients in the asymptotic Bondi 
coordinates are given by ||l7|,[lq| 



e" =z;2(sim?){l + 0(l)}. 


(44) 




(45) 


p^^ = (sin^){l + ^ + 0(l)}, 


(46) 


7 1 
P^v^ - + O(^), 

V V 


(47) 


2m , 1 , 

2/1 = 1 + o(-), 

V V 


(48) 


A- = O(-), 

V 


(49) 




(50) 



where the expansion coefficients a,/3, 7, and m are functions of [u^'d,ip). Then the total energy at the null infinity 
coincides with the Bondi energy C/b(u)j 

lim U{u,v) — f dfl m{u,'d,ip) = ?7b(u), (51) 

where m{u, -d, (p) is the mass aspect of the asymptotically flat radiating spacetime. One can easily show that the 
equation (|9|) becomes 
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^Ub{u) = --^ (fdnv^p^''p'''{d+pac){d+Pbd) < 0, (52) 

which is just the Bondi mass-loss formula. Notice that the negative-definite energy-flux is a bilinear of the traceless 
current j\ defined as 

]\ p'^'d+p,, if, - 0), (53) 
representing the shear degrees of freedom of gravitational radiation. 

III-2. The Misner-Sharp energy 

Let us now consider a spherical ball of radius v filled with a perfect fluid with energy density p(v). The Einstein's 
equations are modified by the presence of the fluid, and the solution is given by 

ds^ = -2dudv - 2h{v)du^ + v^dn'^, (54) 

where 

2hiv) = 1 - (55) 

V 

m{v) ^Airf dv'v'^ p{v'). (56) 
If we choose the subtraction term such that 







then U (w) becomes the Misner-Sharp energy m{v) ||l9| , ^ 

U{v) = m{v). (58) 



III-3. Black holes 

One might be also interested in applying this formalism to black holes, and try to obtain quasi-local energy of the 
event horizon and quasi-local energy-flux incident on the horizon. For this problem, it is appropriate to choose a 
coordinate system adapted to the m-going null geodesies. In a spacetime where the metric is given by 

ds^ = +2dudv - 2hdu^ + e^pab {dy" + A_^'du + A^dv) {dy^ + A^du + A^^dv) , (59) 

the vector field 

is an in-going null vector field generating the null hypersurface u — constant, and the vector field 



(61) 



whose norm is —2h is either timelikc, null, or spacelike, depending on the sign of 2h. Thus, in a region of a spacetime 
where 2h ^ 0, it becomes an out-going null vector field. Therefore, on the event horizon H generated by the out-going 
null vector fields, we must have 2h = 0. If we repeat the previous analysis in the in-going null coordinate system 
(p9|), we obtain another set of quasi- local conservation equations. On the event horizon H, the quasi-local energy 
conservation equation becomes 

^ " i £d'y{le'^p'"'p'\D+pac)iD+pM) - e^iD+af - 2e^«;i?+a}, (62) 

Un:^-^ £d^yie''D+a) + U°i, (63) 
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where n := D_/i|h is the surface gravity on H. This equation is identical to the quasi-local energy conservation 
equation on the stretched horizon, which was studied in detail in ||2l|] . Notice that when the subtraction term is 
chosen zero, the quasi-local energy C/h is non-positive since the area of the event horizon always increases, 

D+a\n > 0. (64) 

However, when the black hole no longer expands so that D+ctIh = 0, then U-a becomes zero. For instance, for 
a Schwarzschild or Kerr black hole, we have C/h = This counter-intuitive aspect is a manifestation of 

the well-known teleological nature of the event horizon. That is, when the event horizon H evolves, its quasi-local 
energy must be negative so as to cancel out the positive in-flux of energy carried by subsequently in-falling matter 
or gravitational radiation, leaving J7h = when the black hole reaches the final stationary state. Details of this 
derivation and discussions of the remaining quasi-local conservation equations on the event horizon will be presented 
elsewhere fl^ . 



IV. DISCUSSIONS 



In summary we derived a set of four quasi-local conservation equations from the Einstein's "constraint" equations. 
In particular, we showed that one of the quasi-local conservation equations reproduces the Bondi mass-loss relation 
in the asymptotic region of asymptotically flat spacetimes, and that the quasi-local energy coincides with the Misner- 
Sharp energy for spherically symmetric fluids in a finite region. We also applied the quasi-local energy conservation 
equation to the horizon of a general dynamic black hole, and found that it reduces to the quasi-local conservation 
equation of Thorne et al. It must be stressed that our quasi-local conservation equations are exact and unique, in that 
they were obtained directly from the Einstein's "constraint" equations through the first-order canonical formalism. 

It seems appropriate to mention that the equation ( pO| ) has a similar structure to the integrated Navier-Stokes 
equation for a viscous fluid [p5|| , 

= - j> dS^ {p6ik + pv^Vk - a[f^ , (65) 

where Pi and a[^. are the total momentum and the viscous term, 

P,^ jdVipvi), (66) 
, (dv,dvk 2 dvi\ dvi 

and T] and ^ are the coefficients of shear and bulk viscosity, respectively. This equation tells us that the rate of the 
net momentum change of a fluid within a given volume is determined by the net momentum-flux across the two- 
surface enclosing the volume. Notice that the Hamiltonian function H in ( p^ ) is at most quadratic in the conjugate 
momenta tt/, and assumes a form of momentum-flux of a viscous fluid. From this point of view, terms quadratic in 
TT/ are responsible for direct momentum transfer, terms linear in tt/ may be regarded as viscosity terms, and terms 
independent of tt/ as pressure terms. This observation allows us to interpret the Hamiltonian function H as the 
gravitational momentum-flux and the two-surface integral 



1 

16^ 



d'y {nh) (68) 



as the quasi-local gravitational momentum associated with iV2. 

The equation (|3l|) is a quasi-local angular momentum conservation equation, since the r.h.s. assumes the canonical 
form 

(69) 

representing the angular momentum- flux associated with the vector field ^ = The proposed quasi- local momen- 

tum and angular momentum conservation equations will be analyzed in detail in a separate paper [|l6|. 
A final remark concerns with quantum gravity. By replacing the conjugate momenta tt/ as 
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in (^), ([3^), and (^J), one can obtain a set of functional Schrodinger equations. For instance, from the equation (|3^), 
one obtains the foUowing equation, 



167r 



1 




du ' 



(71) 



where the subtraction term was chosen zero. These functional Schrodinger equations seem worth exploring in 
situations where quantum gravity effects are expected to be dominant, for example, near black holes. 
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